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Preface 


This document contains the proceedings of the NASA Workshop on Computational 
Structural Mechanics, held at NASA Langley Research Center, November 18-20, 1987. 
The workshop was sponsored jointly by NASA Langley Research Center and NASA 
Lewis Research Center. 

The purpose of the workshop was to allow participants in Langley's and Lewis’ 
Computational Structural Mechanics (CSM) research programs to meet and to 
share research objectives and accomplishments. The intent was to encourage a 
cooperative Langley/Lewis CSM program in which Lewis concentrates on engine 
structures applications, Langley concentrates on airframe and space structures 
applications, and all participants share technology of mutual interest. 

The workshop was organized into the following three sessions: 

I Concurrent Processing Methods and Applications 

II Advanced Methods & Testbed/Simulator Development 

III Computational Dynamics 

Session I dealt with parallel processing methods and languages, new computer 
hardware, and software architecture to exploit parallel computers. 

Session II dealt with the Langley CSM Testbed, the Lewis Engine Structures 
Computational Simulator, and Structural Analysis Technology involving finite 
elements, boundary elements, and probabilistic approaches. 

Session III dealt with advanced methods for structural dynamics. 

The use of trade names or names of manufacturers in this publication does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


W. Jefferson Stroud 



With the exception of a few adjustments made 
primarily for the purpose of uniformity, all pa- 
pers have been published as received. 

— Editor 
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Current 2D/axisymmetric capabilities 


CL 

o 



c 

0) 


E 

o. 

o 

Q) 

> 

CD 

"O 


Q 

co 

h- 

C/) 

UJ 

CD 


C/> 

a> 

(0 

> 

CD 

c 

CD 

a 

CO 


(/) 

LU 

CD 

"O 

a) 

+-» 

o 

o 

o 

<0 


(/) 

C/> 

a) 

o 

o 

a 



CD 

c 

CD 

s 

c 

CD 

'35 

a) 

"U 

c 




o o o o 


764 



I 0 ) 

x 

< -m 


c: 

• o 



m 2 
u 
a) o 

X MH 


o ,r * 

<D J* 

•iH i — \ 
> (0 
0 ) -P 

u 

<D 

<u x 

> -p 


to ^ 
C *0 
0 ) o 
x x 

0 ) -P 
U d) 

a & 

£ 

O < 1 > 
U X 

-p 

(0 


0) o 
> 

■h c/3 

co 
O (0 
-P X 


- a? 

^ x 


H 

4J 

(0 


-p HH 


o 

4h 



> 

• »H 

0) . 

u 

•iH (JJ 

X 

> T3 

fl) »H 

(TJ 

U r—i 


CO 

C 

-p 

•iH 

U (0 
O “'H 


X X 

<1) 

co -p 

rH 


X 

>1 c 

«*H 

rH o 

CO 

(U 

co 

S T5 

o 

a <u 

a u +j 


-P CO 

-P 

X -h 

0 

CD rH 

c 

C CO 

CO 

co O 

• rH 

• iH 


u a 

-P 

<U 0 


•p -p 


765 



FOR A HOMOGENEOUS ELASTIC STRUCTURE 
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The BEM is based on the application of the Reciprocal Work Theorem with (for 
isotropic materials) the Kelvin point load solution and its derived traction solution 
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. INTERIOR DISPLACEMENT EQUATION 
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Primitive Displacements Displacements 
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The computing and modelling characteristics of the BEM are fundamentally different 
from those of the finite element method. For elastic analysis the BEM models both 
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2D/AXISYMMETRIC BEM CAPABILITIES 
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body forces 
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Advanced capabilities are presently being incorporated in the Pratt & Whitney code, 
and will become available in the production version over the next year. 
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BEST 3D is a general purpose BEM structural analysis code developed by Pratt & Whitney 
and the State University of New York at Buffalo under NASA contract NAS3-23697. 
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substructuring for all problem types 
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Major advances have been made during the development of BEST3D, both in the creation 
of new BEM analysis capabilities and in the incorporation of available capabilities 
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of BEST 3D is the provision of a very complete plasticity capability 
d analysis. In BEM analysis of problems with material nonlinearity 
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PERFORATED PLATE 
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ssical problem of a tension loaded perforated plate was analyzed using the two 
BEST3D model shown. Plastic strain distribution was modelled only in the re- 
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Finite Element Models for Benchmark Notch Analysis 
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Boundary Element Models for Benchmark Notch Analysis 
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BENCHMARK NOTCH - FIRST CYCLE RESPONSE 
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NOTCH STRAIN - MICROSTRAIN 



Use of the weighted mesh shown on the previous chart gave good agreement with exp 
imentally measured strains at the notch root, over two and one-quarter load cycle 
Use of the finest finite element model gave equivalent results, in about twice th 
CPU time. 


I 

U • 

0) (D 0) 


803 



RESIDUAL STRESS IN QUENCHED BAR 
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A problem of significant interest in the gas turbine industry is the development of 
residual stress in large parts during forming and heat treatment. BEST3D was used to 
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NONDIflENSIONAl RADIUS - R/I 2S 


The temperature gradient was initially very severe (over 1000° F between the surface 
and center of the specimen), dropping to essentially zero at large times. 
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8EST3D - 6 CELLS 



BEST3D analysis shows good agreement with the experimentally measured stresses, 
for the relatively coarse model used. 
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Goal is to calculate effectiveness of 


CO 


0 

> 
■ mmm 

+-» 

’co 

c 

0 

CO 





c 

CO 0 

0 g 

o Q- 

ii 

4 -» 0 


T 3 

0 

CO 


CO 

’co 

15 

c 

cu 

o 

■ mmm 

c 

o 

£ 

cu 


0 

E 

■ MM 

I— 



> o o 


810 



To calibrate the dynamic capabilities of BEST3D a BEM model was constructed to evalu 
ate the effectiveness of trenches for the vibration isolation of sensitive equipment 
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Passive Isolation Testing 
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VIBRATION EXCITER FOOTINGS 



Extensive test data is available for this problem. The relative locations of the 
exciter, trench and vibration pickups are shown on this chart. 
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ion Factor Contours 
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BEST3D predictions of amplitude reductions and amplifications show good agreement 
with the experimental results. 
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COMPRESSOR DOVETAIL ANALYSIS 
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Pratt & Whitney is presently in the process of applying its interactive 
dimensional BEM system to the determination of stress concentrations in 
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The new BEM system shows excellent agreement with the old system for a variety of ac 
tual engine geometries. 
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Local implementation should be a 2-4 week 
process 
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NOMINAL 10, 0 KSI 
MRX 27, 9 KSI 
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It was possible to provide a 3D comparison between the baseline geometry and 
posed modification rapidly enough to guide the design process. 


I 

I Q • XX rH I 

O co CO U 03 X 

P DO 0) *H 

Qh l/) CO i — I i — I 

O >1 c/3 

(13 H CO 03 g > 

W (0 (U-H O 

■H c ^ 

w nj a) 

>1 p h o> 

•h x a -«h p p 

oj a; c Xi 03 p 

ch ^ td > u 

(0 CX P o< *h • 
g CO 03 4-1 4-4 CO 
U O U O 03 

-rH O 4-1 "O ■ rH 

P O OIW C 

CO CO C *rH f0 -rH 

<13 C/3 CO -rH CO H 

•P 03 -rH +J >i £ .rH 

03 rH CO 3 i-H P X) 

>1 Ql 03 O CO 

03 03 rH g c 4-1 Q« 

C,c (0 o (0+) <0 

•rH -p C U (DO 

P CO P H 

2 +j (0 acn 

OO4C03 C 
P 4-H (0 03 c: Cn -rH 
03 CO -rH C CO 

<D C O H -rH C/3 
1 — I • — 1 -rH ^ C D 03 

U £) H Q 4 O O U 

•rH (0 C C rH O 

(0 h o cn op 

4-1 -rH G C 03 c Oh 

<0 -rH Cj -rH 
^ > * (0 -rH rH 

O 05 CP C P 03 
•P C 2 03 H 

lH >*4 • rH ^ O T5 rH 

03 rH I— I 03 Oh (0 <0 

-rH I — | f — ( * — | p 

•H CO 03 £3 X3 (0 

rH (0 TJ *rH • t 

•rH 03 o CO CO 03 
> g <0 -rH C 4-1 

T3 03 CO -rH O 
CO C 03 4-1 co X3 

03 3 > X5 P C 

•rH O *rH 03 P O 

P P P X) P P *rH 

-rH (0 U x: P 

H c 03 o CP 03 (0 

-rH P 4-1 CO -rH P p 

X) O 4-1 H c -rH -rH 

(OP <D 03 4P O 

Q* 03 C rH 

(0 03 o rH > 03 CX 

U g P H O X 

•rH -rH C 03 

P P C •* C 03 
03 >1 O 03 03 

P 03 -h (0 CO X: 
p'd PH C 03 P 
Q 4 C«rH O 

£ x: 03 o 4 =: o; 

O P P 4-1 >4 U P 

U -rH P PH O -rH 
P > 03 -rH X) (0 D 

03 03 03 tj 1 

a - p X} - au 

2 CO 03 03 O CO P 

CO 03 O4 C p 03 

H O -H 04 P >1 

P P P X) 2 h 

cp apx px 

03 03 P CPU (0 

co g x p p c X) 

03 o p O p o 

P 0)*H COX P P 

CU cn£ o) p co Q 4 


829 




830 


DEVELOPMENT OF AN INTEGRATED BEM FOR HOT 


N89- 29801 


& 

o 

M 

H 

O 

W 

H 

fc 

t—i 

H 

P 

Eh 

0 

P 

H 

03 

1 

Q 

& 

Ph 


v CO 

S3 So 

9 a 

« p 

d 


“'s 

•s « 

<L> -u 

g «3 

•S ^ 

bO m 

H 

> <D 

b £ 

t4-H 

o 0 

-4-9 ^ 

fl .15 
0 ) co 

a s 

H .SS 
3 rt 

Qi 

a> ^ 

p— i ^ 

cd 

49 

C /3 


a 

a 

p 

03 


PRECEDING PAGE BuANK NOT FILMED 



'// y- ? 


! ;? 


o c i 


831 


» Mg ^ Iff J MTBfflONAUl 81AHI 



DEVELOPMENT OF AN INTEGRATED BEM FOR HOT 
FLUID-STRUCTURE INTERACTION 




CO 









H 









55 









W 








W 

2 



H 

s 

1 





0 

w 







►H 

H 

C5 

. HH 




55 

a 



O 

< 

a; 

a. 

D 

cr 

W 

OS 


CO 

w 

w 

X 
H 
1 * 

X 

0 

hH 

HH 

H 

< 

H 






> 

U* 

H 

X 

CO 


r-* 

•—4 


HH 

O 

< 

W 

w 


w 

* 3 ; 

0 

H 

H 

O 

1 

X 

*4 

X 

2 

w 

*—3 

CU 


C 

H 

55 

W 

5E 

CU 

O 

U3 

H 

2 

i-4 

2 : 

55 

O 

§ 

>* 

< 

•n 

CQ 

O 

< 

H 

CO 

X 

O 

Uh 

CU 

55 

HH 

«c 

K 

w 

H 

O 

U 4 

0 

:z; 

< 

►J 

w 

* 

* 

* 

* 

* 

D 



> 






Q 

3= 

W 

w 






O 

W 

05 

Q 






X 

►H 

D 







H 

> 

H 

2 






55 

w 

D 

w 







a: 

U« 

CQ 






* 

* 

* 

* 







>< 

a: 

< 

2: 

2: 

D 

CO 

* 


832 


FUTURE DIRECTION 


>* 


H 

2 

W 


N 

2 


a 

< 

•j 

pa 

H 

fa-i 

W 

hJ 

W 

a 

< 


CO 

M 

K 

E-* 


833 






835 



FUTURE ANALYTICAL REQUIREMENTS 
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VISCOUS FLOW 


The BEM has reached a very high level of j^rformance in the solid mechanics 
where elastic, elastoplastic, viscoplastic, thermoelastic, thermoplastic 
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BOUNDARY INTEGRAL EQUATION 


VIEW GRAPH 7 

One of the biggest disadvantages of BEM is that it cannot be developed if 
the fundamental solution for the set of governing differential equations of 
the problem does not exist. Unfortunately, this is true for the coupled 
set of governing equations for the present problem. 
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FIGURE 1 - Dtivm Cavity Ptobl« Definition 
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FIGURE 5 - Driven Cavity (Three Region, Thirty-six Cell Model) 
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DRIVEN CAVITY (INCOMPRESSIBLE VISCOUS FLOW) 
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HORIZONTAL VELOCITY at X=0.5 
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An outline of the schema in applying probabilistic finite elements is shown on this slide. 
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The starting point for the development of the probabilistic finite element method is the Hu- 
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SECOND MOMENT METHOD OPERATORS FOR 
FORMULATING PROBABILISTIC FINITE ELEMENTS 
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The development of the probabilistic finite element method is based on a second-order moment 
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The governing equations for the system can then be written as follows. The mean value 
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ITEM CALCULATED QUANTITIES EQUATION NUMBERS OF 

NUMBERS INTEGRATIONS 
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The number of equations to be solved is summarized in this slide. The mean values are 
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We will now consider several sample problems which have been solved in order to illustrate 
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The two-degree-of-f reedom problem was solved by three methods: 
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This slide shows the variance in the displacement at node 1 as predicted by the three 
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This example is far more complex and involves strong nonlinearities because the yield stress 
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This shows the mean displacement at node 1, which again has been calculated by the 
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This shows the variance in the displacement at node 1 for this nonlinear example 
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This shows the way displacement output would be presented to a typical user of 
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This slide shows a similar representation for the stress in element no. 1. This can be seen 
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This figure shows a simplified model of a turbine blade which has randomness both in the 
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This slide shows the expected value of the displacement at node a, labeled "expectation 
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Note: Zeroth Order Eqn. is Nonlinear Whereas 

the First and Second Orders are Linear 


To provide some guidelines as to why Monte Carlo simulation cannot be used for problems of 
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PROPOSED WORK 
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This slide summarizes the work which is proposed in the continuation of this grant. The 
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p n [f] _ rf(1,s,t) + (1-r)f(0,s,t) 

P S [F] = sf(rXt) + (1-s)f(r,0,t) 

Pj[fJ = tf(r,s,1) + (1-t)f(r,s,0) 
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FORMULATION OF A TRANSITION ELEMENT 
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A 12-ncde transition finite element 
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By utilizing similar interpolation and blending functions various size finite element meshes are 
obtained from a closed shell structure. 


997 




998 


Iwo (lLmcnsionnl bending finite elements based on the boolean sum are utilized for solution of a 
rectangular plate under uniformly distributed loads. The results obtained are in excellent 
agreement with the theoretical values even with a coarse mesh size. 
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DESCRIPTION OP A BODY IN MOTION 
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Elevated temperature effects, centrifugal forces and time dependent mass effects will be incorpo- 
rated in the so 1 ut i on . 
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RECENT PROGRESS - Out-of-House 
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LATDYN CAPABILITIES 
Finite Element Based 

• General geometry including closed loop topologies 

Rigid or Flexible Beams 

Applied Forces. Displacements. Velocities or Accelerations 

Fortran Based Command Language 

• Permits external user access to internal fortran code 

• User written logic commands provide added control 
over program flow 

Nonlinear Springs and Dampers 

Control Forces 

• Functions of system variables 

• Time delay 

Lock-up, Docking. Impact 
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LATDYN CAPABILITIES 



x 

, 












2? 





73 



4-) 

g 




• 



73 





g 





d) 



•H 

<d 


73 


rH 


• 

d) 


73 



•H 

d) 



U 

CO 

• 

d> 

2 

M 

u 

d> 

CD 

d) 


>1 



0) 



P * 



o 

G 

>1 

u 


oi a) 

N 

N 

> 


p 

£ 


*H 




id 73 


p 


rH 

<d 73 

o X 

•H 

•H 

d) 


-H 

P 


rH 




g 

d) 



CO 

G 

x 

a> 

H 

p 

rH 

U 

rH 


G 

H 

g 

a 



CO 


73 


g 

-H 

£ 

z 

p 

a cd 

fd 

d) 



g 


fd 

a 



d> 

W 73 


X 


73 


fd 


M 

u 

P 

P 



73 

d) 

fd 


<D H 73 

fd 


p 


T5 

to 

-H 

p 


d) 

O 

C 

*» 

£ 

C X 



X X 

G 



•H 

p 

G 

g 

O 

a 


G 

fd 

d) 

>i 

o 

fd 


u 


P 

fd 

fd 

G 


H 

H 

co 

0) 

o 

a> 

CO 

d> 

M 

G 

H 

o 


d) 

d) 



•H 

g 

id 


O 

fd 

d) 

1 — 1 

CO 

g 

d> 

Oi 

fd 

O 

P 


CO 

X X 


Oi 

>H 

£ 



Oix 

X 

CO 

a) 

-H 


X 

a d) 

to 

M 

p 

p 


G 

fd 

o 

M 


rH 

o 

Oi 

o 

fd 

i — i 

X 

** 

o 

g 

g 

d) 

d> 


-H 


-H 

> 

o 

d) 


< 


G 

M 


a) 

o 

to 


o 

O 

u 

a 

CO 

d) 


P 

g 


CO 



p 

fd 

a >i 

X 

d) 

(0 

o 

d) 

G 

g 

P 



P 

o 

G 



G 

X 


oi a) 


X 

c 


cnp 

fd 

’d 

>i 


-rj 

a> 

-H 



CO 

d> 

o 

4-1 

u 

p 

p 

o 

o 

<d 


o 

73 

£ 

P 


g 

p 

OI d) 


a 

a 


O 

a) 

*H 

G 

o 

•H 

X 

rH 

G 


M 

-H 


H 

CO 

O X 


H 

cr 

CO 


a 

G 

d) 


p 

p 

fd 

H 

Ok. 

<P 

O 


d) 

>1 rH 

P 


o 

<u 

p 

CO 

0) 

-rH 

G 

p 

o 


u 

P 

CO 

a d> 


a 

CO 




o 

co 

G 

0) 

P 

O 

G 

G 

G 

d) 

C0 

Oi 


a 




d) 

P 

• 

< X 

*H 

CO 

c 


a 

a> 

G 

o 

C 


G 


(0 


>ip 

CO 

•H 

C 


G 

fd 

CO 

-H 

>1 g 

g 

P 


d) 

d) 

•H 

p 



X 

o 

«H 

g 

o 

CO 

CO 

M 

fd 

<d 

G 

o 

0) 


CO 

o*x 

M 

u 

>1 


(U 


fd 

H 

•H 

fd 


P 

rH 

H 

(d 

o 

1 — 1 

G 

d) 


p 

a 

fd 



G 

CO 

p 

<D 

P 


fd 

CO 

O 

P 



<u 

o 

a 

H 


CO 

X 

g 


d) 

c 


a g 

G 


C 


CO 

c 

d) 


•H 

fd 

O 

G 


o 



X 

o 

u 


O 

S 

G 

O 

c 


-H 

X 

d) 

P 

X 

P 




H 


P 

-H 

o 

o 

d) 

o 

O 

o 

-H 

a) 


p 

p 

fd 

CO 



CO 

p 

d> 



P 


p 

X 

G 




g 

CO 


•rH 

>H 


d) 

>1 

g 

G 

CO 


CO 

o 

d) 


d) 

73 

g 

0) 

o 

fd 

p 

G 

O 

rH 

o 

P 

rd 

d) 

G 


X4 G 

G 




c 

V 

CO 

CO 


o 

*H 

a 

fd 

G 

•H 

d) 

o 



-H 

G 

H 


p 

g 

fd 

p 

3 


CO 


P 

u 

X 

fd 

M 

X 

a 1 

d) 


X 

P 

P 

• 

o 

X 


CO 


CO 

P 

G 


d> 

o 


fd 


d> X 


CO 



o 


P CM 

>1 

Jh 

•H 

G 

o 

P 

p 

rH 

o 

rH 

d) 

CO 



G 

d) 

G 

p 

H 

•H 


CO 

O 


0) 

-H 

CO 

G 

O'* r — 1 

G 

rH 

X 



O 

M 

d) 

d) 

d) 

M 

d) 


4-1 

a) 

g 

p 

o 

-H 


i — i 

a 

X 

G 



•H 

fd 

P 


73 

O 

o 

<D 


M 

0) 

o 

g 


M 

fd 

O 

•H 

CO 



P 


P 



O* G 

JH 

73 

0) 

rH 

g 


d) 

O 


a 

X 


• 


fd X 

*H 

Du 

O 

<H 

d) 

d) 

a) 

X 

d) 


c 

ap 

CO 


d> 

fd 

CO 


rH 

o 

H 

H 


fd 

M 

X 

73 

p 


d> 

-H 

fd 


P 

d> 

rH 


c 


d> 

-H 

? 


H 


d) 

s 

C 


d> 

N 


X 

d) 

-H 

73 

P 

G 

o 


M 

X 



-H 

G 

p 


d) 

«k. 

P 

•H 

CO 

CO 

73 


•H 


O 

-H 



S 

VH 

• 

d) 

fd 

(U 

G 

P 

CO 

*H 

H 

< 


fd 

CO 

3s 

H 


P 


73 


d> a x 


u 

O 

a 

*H 

G 

d> 


73 

g 

■H 


O 

73 

fd 


d) 

CO 

CO 

O 

P 

73 


•H 

*H 


* rH 

P 


G 



CO 


d> 

g 


G 

d) 

G 

• 


C 

C 

P 


P 

P 

O 

• 

fd 

CO 

X 

p 

73 

X 

M 


•H 

o 



73 

fd 

*H 

fd 

CO 

fd 


fd 

CO 


•H 

o 

-rH 

•H 

-H 

o 


P 

M 

w 


G 


73 

G 

*H 

X! 

P 

u 

p 

co 


fd 


Oi 

H 

p 


d) 

O 

d> 

• 

fd 


P 


P 

O 

fd 

G 

P 

G 

o 

>1 

-H 

O 

d) 

• 

73 

P 


Q 


73 

a> 

•H 

4-> 



X 

0) 

C 

O 

u 

X 

H 

CO 

73 

CO 



9 

s 

73 

d> 

CO 

CO 

M 

u 

X 

o 

g 

•H 

-H 

a 



d) 


G 

H 

< — i 

g 


d) 

P 

CO 



c 

CO 


d) 

O 

CO 

a 

o 

O 

73 

d) 

O 

d) 

o 


• 

P 

P 

G 



o 

d> 

o 

rH 

a 

•H 

fd 

p 

P 


01 

-H 

CO 

M 

d) 


G 

•H 

O 

d) 

• 

*H 

g 

p 

d) 


> 




tat 

M 

P 

G 

P X 


O 

g 

CO 

X 

g 

P 




X 

O 

P 

73 

73 

g 

fd 

fd 


G 

p 

CO 

d) 

M 

♦H 

P 

m 

fd 

fd 

CO 

d> 

G 

u 

G 

<V 

d) 

d) 

rH 

H 

P 

O 

M 

fd 


d) 73 

o 

M 

g 


dJ 

P 

O 

Or d) 

P 

P 

p 


d) 

O 

O 

G 


d) 

a 



O* 

H 

X 

a 

•H 



g 

CO 

•H 

CO 

73 

H 



|X4 X 



CO 

>1 

O 

O 

P 

<d 

G 

P 


<D 

d> 

g 

>1 c 

d) 

P 

w 


o 

d> 

<u 

fd 

G 

M 

P 

-H 

X 

-H 

G 

M 

H 

p 

•H 

CO 

fd 

U 

G 

CO 


G 

M 

M 


fd 

a 

a) 

5 

CO 

P 

d) 

d) 

d) 

p 

rH 



O 

a 

u 

■ 

CO 

fd 

fd 

V 


73 




g 

> 


fd 


d> 

CO 

fd 

G 

d) 

CO 

■ — ■ 




M 

*H 


73 

d> 

O 

d) 

d) 

d) 


CO 

p 

G 


•H 

a <u 


CO 

co 

CO 

O 

fd 

<u 

d> 73 

P 

H 

S 

p 

CO 

•H 

fd 

O 

H 


g 

-H 

s. 

73 

P 

g 


M 

g 

I - 1 

o 

G 

d) 

O 

-H 

fd 


G 

-H 

O 

<D 

fd 

P 

CO 

c 

G 

d) 

P 

G 

o 

<v 

g 

-H 


X 

G 


<0 

-H 

P 


X 

73 

*H 

p 

fd 

♦H 

p 

o 

P 

CO 

73 



d) 


-H 


P 

73 

fd 

CO 

p 


rH 

G 

g 

fd 

co 

fd 

O 


o 

G 

73 X 


P 

CO 

G 

H 

P 

d> 



-H 

9 

£ 


>i a 

G 

CO 

g 

o 

0) 

P 

* 


d) 

d) 

o 

o 

-H 

CO 

CO 

X 

g 

O 

p 

CO 

g 

M 

d) 


a 

N 


P 

fd 

-H 

g 

O 

H 

P 

5 

O* fd 

d) 

O 

CO 


•H 

4-> 

o 

CO 

G 

•H 

p 

G 


O* 

d) 

O 


-H 

O 

G 

ap 


G 

o 


CO 

Oi 

-H 


P 

o 

d) 

p 

O 

rH 


73 

O 

H 

-H 

fd 

fd X 

O 

-H 



H 


73 

d> 


g 

o 

rH 

d) 

73 

d) 

O 

rH 


o 

P 

O 

o 

£ 

OI 

P 

0) 

d> 

c 

U 

G 

d) 


O 


(D PH 

fd 

a 


CO 

♦H 


fd 

G 

G 73 

M 

d> 

O 

O 

rH 

d> 

a 

P 

P 

-H 

d) 


CO 

u 


X 

Oi G 

•H 

d) 

g 

G 

a 

CO 

-H 

CD 

Oi o 

o 

O 

g 

> 

d> 


d> 

o 


G 


g 

G 

P 

d> 


P 


fd 

p 


d) 

•H 


Oi 

M 

X 

•H 


-H 

a 

O 

0) 


O 73 73 

O 

d ) 

p 


>1 > 

rH 

V 

fd 

fd 

p 

Oi n 

>i 


O 

rH 

g 

G 



g X 

G 

a 

H 

G 

G 

CO 

G 

d> 

o 

o 

d> 

U 

C 73 

a) 

a) 

H 

p 

d) 


P 

fd 

o 

fd 

O 

G 

p 

Oi G 


rH 

> 

fd 

H 


p 

P 

o 


d) 


> 

o 

M 

O 


G 

G 

-H 

p 


o 

> 



a) 

CO 

CO 

c 

<d 

X 

O 73 

rH 

X 


O 

<D 

fd 

r — 1 

o 




CO 

a 

P 

>i 




p 

P 

fd 


*H 

fd 

Oi g 

• — i 

C 


<1 

rH 


d) 

G 

*H 

CO 

<U 

C0 

CO 




73 

rH 


U 

CD 


o 

p 

Pu 

o 

g 

OI 

i 

C 


X 

d> 

d> 

Xl rH 

P 

<D 


p 

d) 

O 

73 

G 

CO 

Eh 

u 

H 

G 

•H 

>iEH 

o 

•rH 

fd 

(d 

G 

CO 

d) 

O 73 

fd 

d> 


o 

oc: 

p 

P 

fd 

o 

4-1 

73 


73 73 

g 

G 

id 

O 

X 


G 

1 — 1 

CO 

P 

X 

O 

G 


X 

o 


O 



o 


M 

o 

«H 

p 

d> 

G 

Or 

fd 

O 


(X4 

O 

«» 

o 

t — 1 

Id X 

■ 

g X 

CO 

d) 

•H 

O 


CO 


CO 

X 


fd 


O 

G 




1 

G 

fd 


<P 

P 

P 


V. 

D 

O 

•H 


G 


rH 


O 

p 

g 

CO 

•H 

O 

u 

p 

P 

G 

*H 

01 

>1 


P 73 


O 73 rH 

P 

-H 

G 

H 

-H 

P 

-H 

O' G 

fd 

•H 

G 

G 

rH 





-H 

c 

id 

O 

P 

-H 

o 


rH 

P 

O 

d) 

G 


Oi 

-H 

P 

• 

C0 


K 

P 

fd 



-H 

(d p 

55 

G 

O 

M 

c 

-H 

G 

-H 

73 

G 

CO 

M 

CO 

Eh 

-H 


P 

d) 

73 

u 

M 


g 

g 

a o t) 

O 

C0 

G 

d) 

d) 

dJ 

d) 


G 

>i o 

d) 

73 

P 

d) 

a 




a, 

M 

•H 


rH 

CO 

CO 

Q 

0 

O 

•H 

fd 


M 

fd 

to 

a 

Eh 

a) 

d> 

d> 

g 

O 

P 

a> 

O 

d> 

CO 

£ 

M 

Cli 

P 

rH 

d) 

Oi 


c 


< 

x 

X 

X 

8 

O 

O 

G 

G 

H 

fd 

d) 

O 


d> 

d) 

CO 

d> 

c 

0 

o 

3 

p 

P 

EH 

o 

O 

g 

o 

•H 

a. 

g 

g 

P 

< 73 73 

G 73 

M 

u 

p 


1021 





p a) 

p p 

-h p w c 

_ c g > 

CO 73 O CD t0 

<d g p p 

4-J (0 ro w 

(0 M >,4J P 

o co to co g a) 

-H 73 O O CO 

C C to M JC D 

G m cd p 

g g t3 p -h a) 

g g a g s tn 

o o to a co 

O O g CD P 

CO O TJ d) 

p a) j* o o > 

CD P O O CO 

CO p O <D 

G XXH <D 

to p m p 

CD P c -P 

P 0) • co P 

H P C0*H CD P 

CX CD > P O 

P **H . G P 

■ CD P g *H 

g 4J *H (0 0) 

(0 G rH M CD g 

P -H *H tnp o 




X) o 

p 

C0 


o 

P (0 P 


P 


p 

o a a o 

CD 


a co to 

p 

p 



CO O CD 




£ 

CD P 

CO 

G 


* 

O g P 

CO 

O 


Q 

O to 

CD 



Eh 

p p p 

o 

>i 


< 

a o 

O 

P 


3 

CD o 

fO 

CD 



p p p 


> 

o 

0) 

a a p 

>i 


& 

p 

<0 

CO 

CD 

w 

p H a to 

P 

s 


fO 

CD 

to 

2 

p 

♦H CD 



Cu 

o 

• p g 

p 

>i 

kl 


CD G O 

CD «H 



tMl) O 

co p 

£ 

p 

t0 P CD 

G 

to 

>H 

0 

G OP 


*H 

Q 

s 

a 

<D 

P 

EH 

Q 

G OP 

(0 

< 

g 

to P P 

P 

> 

Hi 

(0 

H O 


G 


P rH <0 -H 

'W'OHH CD 
g h a) t 3 p 
CD <0 CD P *H O 

p g p co > *h 

p g to ox: 

O CD P S 

co o <o p a 


<D P CO 

-P <D co co CD 

a* X rH 4-> *H G 

P -H rH -H XJ -H 

P H *H g Eh P 

(0 I P P G 

G £ CD O 

h p a • m 

h p 2 o p P 

-H Ej *H g CD 3 

a; p co m w 

>iO SPG 

H tu tJVH c 

H CD O H 0) 

fO (0 73 M P 

o o ax* p 


•h p o 
p & & 

fO G £ >H 


G -H 
2 cd P 


g o Q P 

CD p a H CD p 

P P H *H CD 

O P H -H co 

co o a g 
g ^ 0) g 
CD (0 P O CD 

P P CO p O W 

P tn CD O 

Cr» o p p £ a 

•H P to O ft g 

pap k o 

CD 0 Eh O 

(D (D C P flJ 

P P CD O O O 

E-» P CJ Cn P 


1023 



GENERIC HINGE BODY WITH MEMBER CONNECTIONS 

/ REVOLUTE CONNECTION 



1024 


Define: 

[r], Fixed Orientation of Hinge Connection Relative to Hinge Body 
0, Time Varying Rotation About a Hinge Line 



GENERIC HINGE BODY WITH MEMBER CONNECTIONS 
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It is .important to note that extraction of element deformations only requires defining the convected 
axis joining the two element ends and not any axes orthogonal to it. 


Equations of Motion and Their Numerical Integration 
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EQUATIONS OF MOTION AND THEIR NUMERICAL INTEGRATION 
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COMPARISON OF SOLUTIONS FOR DEPLOYMENT OF FLEXIBLE BOOM 
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ERRONEOUS DESTABILIZING CORIOLIS FORCE ON ROTATING BOOM DUE TO 
COUPLING OF LINEAR STRUCTURAL AND NONLINEAR KINEMATICS MODELING 
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SLIDE 2 


This mesh illustrates why mixed time integration is crucial in the application of explicit 
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SLIDE 3 


This slide shows the fundamental equations involved in an explicit time integration of 
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In order to examine how a multi-time partition is constructed, it is worthwhile to consider 
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SLIDE 5 



When we consider this flow of information in a one-dimensional problem with a time-step 
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numerical experiments indicated that the schemes we have used are stable whenever the elements in 
each sub-domain satisfy a local Courant condition. 
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SLIDE 6 



The computer program has been implemented on an Alliant FX/8, which is an MIMD (multiple 
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SLIDE 7 


This slide shows the flow of computations for a structural dynamics program with mixed time 
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This is the first example we studied. It is an impulsively loaded cylindrical shell. The 
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Cylindrical Panel Timing Study (Seconds of CPU) 
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This gives various timings for the previously described problems. The rows labeled 
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Timings for Containment Vessel Problem 
with and without Mixed Time Integration (in sec CPU) 
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This problem, which is labeled a containment problem, contains elements of various sizes, so 
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REPROGRAMMING IS NECESSARY TO TAKE ADVANTAGE OF 
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INCORPORATION IN TEST BED 


We list our conclusions and suggestions for future work. 
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. Belytschko and N. Gilbertsen, “Concurrent and Vectorized Mixed Time, Explicit Nonlinear Structural 
Dynamics Algorithms” 

- - Paper to be presented at Symposium on Parallel Computations and Their Impact on Mechanics, 
1987 ASME Winter Annual Meeting, Boston, Massachusetts, December 13-18, 1987. 
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‘CUT AND PASTE’ ALGORITHMS 


MODEL STRUCTURE 




THE STRUCTURE VIEWED AS A COLLECTION OF 
DISCONNECTED SUBSTRUCTURES 
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A ‘CUT AND PASTE’ ALGORITHM 


• Predictor phase: 

d n +i = d n + A tv n + (1/2 — /?)A t a n 
v n+ i = v n + (1 - 7) Ata n 

• Equation solving phase: 

&n+ 1 — ^ 

for 3 — 1, NS do 

K+x = — (M s + 
a n +l ^-n+l d - hd ^-n +1 
a n +i M a n+ i 

• Corrector phase: 

d n+ i = d n+ i + f 3 At a n +i 

v n+ i = v n +i + 7 Ata n +i 


fMJ| ***** 
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INTERPROCESSOR COMMUNICATIONS 



REDUCED SUBSTRUCTURES SHOWING THE 
COMMUNICATION DUE TO SHARED DEGREES 
OF FREEDOM . 
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c 


OVERVIEW OF GENERAL PROPERTIES 


• Parameters: 

n = Number of dof in structure. 
s = Number of element groups. 
p = Number of processors. 
i — Number of interface dof. 


• General properties: 

i) Newmark’s method is obtained for 5 = 1. 

ii) Unconditional stability for all 5 and 7 > 

C 1/2, p > 7/2. 

iii) Full concurrency on a p-processor machine 
(p < n) except for 0(i) operation (mass- 
averaging). 

iv) For given accuracy and njs 00 , 


SPEED -UP = 


0(p>/s), ( 2 D) 
0(ps ), (3D) 


MOT FILMED 


toL. 
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General properties. Note two-parameter dependence of speed-up esti- 
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Figure 3. Discretization and partition of the bar problem. 


One-dimensional test for assessing communication efficiencies. 
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Figure 4. Discretization and partition of the plane stress problem on 
32-processor computer. 


Two-dimensional test for assessing communication efficiencies. 


1079 



N 

Computation 
time (milisec.) 

Communication 
time (milisec.) 

efficiency 

% 

Rate 

Mflop 

64 

1.47 

0.46 

76.32 

0.233 

96 

2.20 

0.46 

82.86 

0.253 

160 

3.67 

0.46 

88.96 

0.271 

288 

6.61 

0.46 

93.55 

0.285 

544 

12.49 

0.46 

96.48 

0.294 


Table 5. Performance of the Bar problem on the 32 Processor Hypercube. 


No. of 
Elements 

Computation 
time (milisec.) 

Communication 
time (milisec.) 

efficiency 

% 

Rate 

Mflop 

128 

25.2 

7.8 

76.3 

0.232 

288 

57.1 

8.6 

87.0 

0.265 

800 

161.3 

10.6 

93.9 

0.286 

2592 

530.9 

14.6 

i 97.3 

0.297 

9248 

i 

1915.2 

! 

22.6 

i 

i 

98.8 

1 

! 

0.301 


















Communication efficiencies for 1-D and 2-D test cases. Note efficiencies 
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C COMPUTATIONAL EFFICIENCY 

• COST « \nb 2 + 2 nb , (6 = semi-bandwidth) 

• Square mesh, l 2 elements: 


GLOBAL ~ + 2) 2 (Z + l) 2 + 2(1 + 2 )(Z + l) 2 

• Partitioned mesh, s = m 2 subdomains: 


c 




PARTITIONED 
m J \m 


l 


+ 1 



• Equation solving speed-up (n/s — > oo): 


SPEED - UP(2D) = 
SPEED - UP(3D) = 


GLOBAL 

TartitJoWed 

GLOBAL 

paMttIoned 


0 ( 8 ) 

0(s 4/3 ) 
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2D CASE (1024 ELEMENTS) 


C 



dn-CD3dS 
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0.0 200 • 400 600 800 1000 1200 

NUMBER OF SUBDOMAINS 





Estimated equation solution speed-up for one application of the algo- 
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ACCURACY ANALYSIS 
• Algorithmic phase errors, ID case: 

C 



• Maximum celerity of computed waves 



= AL/At 


A L = subdomain size 

• For accurate results, need to take 



Cmax > c > OR At < AL/c 
wave celerity 
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Accuracy requirements derived from an analysis of phase errors in one 
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C ACCURACY REQUIREMENTS 

• Square mesh, s = m 2 subdomains: 

A L = L/mtt 0(1/Vs) 

At < AL/c - L/mc « 0(1/ \fs) 

• Net speed-up (p = 1): 

SPEED - UP(2D) « 0(s ) x 0(1 / a/s) = 0(/i) 

• Cubic mesh, s = m 3 subdomains: 

c 

A£ = L/m « Oil/s 1 ^) 

At < AL/c = L/mc « 0(l/s 1,/3 ) 

• Net speed-up (p = 1): 

SPEED - UP(3D) « 0(s 4/3 ; ) x 0(l/s 1/3 ) = 0(a) 


I 
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Estimated equation solution speed-ups for a square mesh in large scale 
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NUMERICAL TESTS 


Square membrane, simply supported, subjected 
to uniform initial velocity. 

Finite deflection FE formulation. Triangular el- 
ements: 


T A 2 


2 A 0 


T = tension 

Ao = initial area of triangle 
A = deformed area. 


Quadrilateral elements: 


Parameters: L = 2, T = 1, p = 1. 
Error measure: 


ERROR — 


X 


dt 

Wexact{t) | ^2* 


1/2 




< 

2 

o 

M 

H 

2 

W 


H 

2 


2 

< 


f- 

3 


u: 

a 

< 

cu 

on 


1093 




1 SUBDOMAIN, h = 0.05 







UNIFORM IMPACT ((3 = 0.25, 7 = 0.5) 


( 



UNIFORM IMPACT (0 = 0.25, 7 = 0.5) 
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ACCURACY REQUIREMENTS 
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NUMBER OF SUBDOMAINS 




Actual vs. estimated time step requirements as a function of number of 
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1024 ELEMENT CASE 
NSUB | Secs. Speed-up Theory 
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Actual vs. estimated speed-ups for square membrane problem on a 
single processor. Timings correspond to the equation-solving phase only. 
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COMPUTATIONAL METHODS AND SOFTWARE SYSTEMS 

FOR 

DYNAMICS AND CONTROL OF LARGE SPACE STRUCTURES 


N89- 29809 


5 / 0 - 3 ? 


1 / A 

^ v O i 


O) 

x 

o 

m 


■3 | 

j=i i; 
H « 

£ « 
■« w 
Is £ 

J=! ^ 

o -O 

<3 

S' 3 
;& .2 

O 

< d 
d ^ 


« o> 
2 o 
<u o 
00 


a 

S u 

0^0 
v 'y n 

~ QJ ^ 
J 3 ^ - 

2 a s 

c *< *g 

n t) o 

u g cq 

ra ^ 
O 


»- 

(X, 


o 

6 Q 


no 

§ 

8 

•«-» 

u 

3 

M 

CO 


id 


o 
cx 

.CO . 

Q Sod 


bC 

•g'S 

a s 

a> — 
C o 

•ao 

C «*- 

w ° 
*s £* 
& 2 
§f £ 


<V 

a 

(V 

u 


u 

< 

a: 

h- 

CD 

< 




Ih 

<y 

a 

V 

O 

*d 

o 

a) 

8 

£ 

bO 

3 

CO 

£ 

a 

.2 

*to 

*> 

5 

Cfl 

U 

fi 

C 

>> 

Q 

f "0 


8 

M 

+J 

u 

2 

-4-> 

CO 

a 

o 

* 

0) 

3 

•2 

o 

6 

a> 

■s 

3 


A 


o 

u 

u 

a» 

IE 

3 

O 

OQ 

o 

"O 

ns 

k_ 

_o 

O 

u 

o 

>* 


a> 

w > 

’c 

3 

jd 

o 

k . 

4 m * 

c 

o 

u 

T 3 

C 

P 5 

L/> 

< L > 


U 

3 


CO 

a> 

u 

ns 

a. 

to 


c 

<u 

U 


PRECEDING PAGE BLAt’vK NOT FILMED 


1105 


FA fift /I bJl I NTENJIDNALU BLANK 



RESEARCH THRUST 


<D 

s g* 

9 ^ 


M «*■* 

a § 
-2 O 

> T) 

« C? 
Q rt 


o: Q 


1106 


Center for Structures and Controls. I niversity of Colorado. Boulder, CO 


CURRENT RESEARCH THRUSTS 


"1 


3 


a J2 


o 


te 

o 

■3 

<3 

a> 

w 

0) 


a 

<y 

w 

0) 


0) 

45 


0) 

cx 

o 

o 

m 

<v 

45 

• 4 J 

nd 

a> 

C 

C£3 

<D 

nd 


M 

u 

<D 

% 

jet 

I 




9 


$ O ^ 

3 S g 

a> w 9 

•a " § 

>s-* 

• —t 


£ 

<y 

1 

I 

v 

43 


0 

1 

*d 

i 


o 

8 

o 

O 

a< 

a> 

nd 

§ 

4: 

u 

M 

td 

a) 

GO 

<D 


(0 

o 

"*■» 

cS 


CO 

o 


o> 

m 

s 


a> 

O 

4^ 

V 

fH 

0) 

to 

2 

u 

45 

CO 


O 

u 

u 

a5 


(-4 

•3 5 

O J3 

« .P 


g 53 


45 45 


a 


a . 
-a ° 


* I .9 

g § 


’Td 

o 

42 


a 

«< O 


- 


” J bb 
•3 ^ jg 

5« 

5 
o 

o 


8 » 
’5? .5 

42 rr~4 


5 
(5 « 
rt 42 


6 v 


x> 

<lt 

(X 

rt . 

O HH 


* o 3 
5 g no 

«3ft 
O 0.0 
S o % 
55 'S > 

■« -S 
c « "O 

t -B 2 

h m O 

2 i Si 

§ J 8 

^2 
0) ^ CX 

«5 *M 
>V 45 O 


a ® p 

, w (X 

< •* a 
• 8 o 

CO +* u 


a 

a> 


i 

•8 

45 


o 


Pi 

o 

'd 

a 

nJ 


ed 

4-> 

c5 
a» 

6 
JJ 

a. ■ 

6 o 

M <U 

4 3 


w 

o 

Pi 

O 


1107 


Center for Structures and Controls. University of Colorado. Boulder. CO 



I 


< 

c 

*— I 

H 

2 

W 

H 

2 




2 : 

< 

cc 

H 

fcJ 

"1 


o 

< 


LO 


1108 



OVERALL RESEARCH APPROACH 
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C-NICE: BACKGROUND 



Center for Space Structures and Controls, Univ. of Colorado, Boulder CO 


An applicationa-softwarc architecture may be visualized as a layer of software that helps the development of applications by 
“cushioning” the interface of the programs with the operating system (see Figure). 
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NUMERICALLY STABLE AS t- 



Numerical Integration 

The numerical integration method is based on a Nordsieck-like method. The displacement, 
ocity, and acceleration are defined at an initial time. A modified Taylor series is used to 
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nlinear Equations of Motion 
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Local Structural Condensation 
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F 2L + X L F 2L-1 + Y L F 2L+1 
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MULTI GRID METHOD 
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SPECIAL ASPECTS 
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Gear - gear interactions cause a wide band width. 


Special Aspects 

Special aspects complicate the dynamic analysis of transmissions. Gyroscopic and gear 
forces cause nonlinear lateral and torsional coupling. Gear tooth passing frequencies are 
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Consider gear tooth interaction. For any gear set the line of force does not pass through 
the gear centers. In the case of spur gears any perturbation of the radial force will result 
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Gear Tooth Model 

Consider the gear tooth interactions of a spur gear set. The contact point varies as the 
angle of the gear set varies. Machining errors cause the contact point to move. High torque 
can cause the teeth to bend. The number of teeth in contact varies as the toraue varies. 
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TYPICAL TRANSMISSION 







Typical Transmission 

In a typical transmission there are many gear sets. Each of these gear sets causes one lo- 
cation on the structure to interact with another point on the structure. Thus, far interactions 
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Transmission Dynamic Analysis 

Complete transmission dynamic analyses are rare in the open literature. David and 
Mitchell (1986) have used a modal balance technique. The problem with modal techniques 
is that the nonlinearities cause the set of modes not to.be closed. This results in side bands 
around the tooth passing frequency. Whenever a solution is found, it is not known whether 
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Application of Multi- Grid Techniques 
To Structural Analysis 
Using a Parallel Transputer Array 
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TRANSPUTER USED IN GRAPHICAL DISPLAY, G. K. ELLIS (ICOMP) 
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The approach is to use an “OCCAM” parallel programming language (specifically de- 
igned for the transputer) and the test bed (with electronically reconfigurable architectures) 
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RECONFIGURABLE PARALLEL ARCHITECTURE WHICH CAN DRAMATICALLY REDUCE THE COST AND 



Potential Impact 

The multi-grid method, although used for years in fluid dynamics, now offers a new 
approach to nonlinear structural dynamics. The computing time does not depend on the 
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